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Abstract 

A recent development in studies of random non-Hermitian quantum systems is re- 
viewed. Derealization was found to occur under a sufficiently large constant imag- 
inary vector potential even in one and two dimensions. The phenomenon has a 
physical realization as flux-line depinning in type-II superconductors. Relations be- 
tween the derealization transition and the complex energy spectrum of the non- 
Hermitian systems are described. Analytical and numerical results obtained for a 
non-Hermitian Anderson model are shown. 
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1 Introduction 

The main purpose of the present paper is to review a new development in 
studies of non-Hermite operators. The review is mostly based on the work by 
the present author in collaboration with Nelson [1,2] and the works following 
it [3-14]. In addition, I report a new numerical result for a non-Hermitian 
ladder system. 

Non-Hermite operators appear frequently in dynamical systems as Liouville 
operators. Although less frequently, they also appear in the context of quan- 
tum mechanics as Hamiltonian operators in the Schrodinger equation. A well- 
known example is the optical potential, which is a complex scalar potential 
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that effectively describes multiple scattering and absorption. A non-Hermite 
Hamiltonian can also emerge when one relates a (d + l)-dimensional clas- 
sical statistical system to a rf-dimensional quantum system by path-integral 
scheme or the Suzuki- Trotter transformation [15]. As a classic example, Mc- 
Coy and Wu [16] showed that equilibrium classical statistical mechanics of 
a two-dimensional asymmetric vertex model can be described by imaginary- 
time quantum dynamics of a non-Hermitian XXZ spin chain. In this case, 
the non-Hermiticity of the spin chain is originated in an external field that 
generates a diagonal flow of edge spins of the vertex model. 

The new development reviewed here resulted from introduction of quenched 
randomness into non-Hermitian quantum systems. A derealization phenomenon 
was found for an especially simple class of random non-Hermite Hamiltoni- 
ans even in one and two dimensions. The Hamiltonian contains a constant 
imaginary vector potential and a real random scalar potential. As the imag- 
inary vector potential increases, all of originally localized eigenfunctions get 
delocalized one by one. One of the remarkable features is that a complex 
eigenvalue of the Hamiltonian indicates the derealization of the correspond- 
ing eigenfunction. Thus we can study the derealization phenomenon simply 
by investigating the energy spectrum of the non-Hermite Hamiltonian. 

The derealization phenomenon has a physical realization as flux-line depin- 
ning in type-II superconductors. Basic correspondence between the dereal- 
ization and the depinning is described in the next section. In Section 3, I 
discuss the complex-energy spectrum of the random non-Hermitian system, 
using numerical data for a non-Hermitian Anderson model. A new result for a 
non-Hermitian ladder system is also reported. Section 4 presents an interesting 
application of Mott's variable-range hopping to the non-Hermitian quantum 
mechanics. 

Non-Hermite matrices with randomness have recently attracted much atten- 
tion from other viewpoints as well. The spectrum of the Fokker-Planck op- 
erator with a random velocity field has been studied in Refs. [17,18]. Non- 
Hermitian random matrix theory has seen great progress recently [5,13,19-26]. 
The present review, however, does not cover these topics. 



2 Non-Hermitian quantum mechanics and flux lines in supercon- 
ductors 
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2.1 Random non-Hermite Hamiltonian and an elastic string in a random 
washboard potential 



A typical Hamiltonian which I treat here is 

Ms M + V (f), 



where p = (h/i)d/dx is the momentum operator, g is a constant real vector 
referred to as a non-Hermitian field, and V is a random potential. A lattice 
version of the above Hamiltonian is given by a non-Hermitian Anderson model 
on a hypercubic lattice, 
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where t is the hopping amplitude, the vectors {e u } are the unit lattice vectors, 
and V$ is an on-site random potential following a probability distribution 
P(Vg). Periodic boundary conditions are applied to wave functions in both 
cases. Note that the non-Hermitian field g plays a role of an imaginary vector 
potential. 

In the Hermitian case g — 0, Eqs. (1) and (2) are reduced to the standard 
Hamiltonians for the Anderson localization. It is widely accepted for g = 
that all eigenfunctions are localized in one and two dimensions. The present 
author and Nelson recently found [1,2] that all of the localized eigenfunctions 
get delocalized one by one as the non-Hermitian field g is increased and that 
appearance of complex eigenvalues indicates the derealization transition. 

As is exemplified in the study of McCoy and Wu [16], a non-Hermite Hamil- 
tonian can have physical relevance when it is mapped to a classical statistical- 
mechanical system with path-integral mapping. By identifying the imaginary- 
time-evolution operator e~ Arn of the Hamiltonian (1) as the transfer matrix 
of a classical system, we can transform [27,28] the matrix element of the time- 
evolution operator between the initial and final vectors, 

-L T H/h 



into the partition function of an elastic string in a (d + l)-dimensional space, 

Z = J Vxe- E *W T W h (4) 
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Fig. 1. An elastic string in a random washboard potential is subject to thermal 
fluctuations. The energy of the string is given by Eq. (5). Through path- integral 
mapping, statistical mechanics of the system is equivalent to quantum dynamics of 
the Hamiltonian (1). 

with the energy of the string given by 



E cl [x(r)] ee / 



dr 




-|.| + V(x-) 



(5) 



The first term of the energy describes the elasticity of the string. (See Ap- 
pendix for details of the above transformation.) 

The energy (5) of the elastic string is re-interpretation of the imaginary-time 
action of the quantum particle. The imaginary-time axis r of the quantum 
system is identified as an additional spatial axis of the classical system. Hence 
the imaginary time L T becomes the system size of the classical system in the r 
direction. The world line of the quantum particle, x(t), is identified as a spatial 
configuration of the classical elastic string subject to thermal fluctuations. The 
temperature of the classical system is given by the Planck parameter h. The 
integral / Vx denotes the summation over all possible configurations of the 
world line, or the elastic string. 

Note that the random potential does not depend on r. Hence the elastic string 
is put on a "random washboard" potential (Fig. 1). The potential tries to 
trap the elastic string in particularly deep valleys and thereby to align the 
string along the r direction. On the other hand, the second term of the en- 
ergy (5), which comes from the non-Hermitian field g, tries to tilt the string 
away from the r direction as is explained below. (In fact, if the random po- 
tential is not present, the energy is optimized when the string is tilted by 
the angle tan -1 ((?/m).) The competition between the above two effects results 
in a pinning-depinning transition of the string from the random washboard 
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potential. This in turn indicates a localization-delocalization transition of the 
quantum particle subject to the non- Hermit ian field and the random potential. 

The reason why the non-Hermitian field g in the quantum system has the effect 
of tilting the elastic string may be explained in the following way. A vector 

— * 

potential A generally induces a current in a quantum system because of the 
coupling p- A in the Hamiltonian. The current would be expressed by the world 
lines of quantum particles running diagonally in the real-time-space. Since the 
present non-Hermitian field plays a role of an imaginary vector potential, it 
induces an imaginary current. The imaginary current is expressed by tilted 
world lines in the imaginary-time-space. Hence the field g tries to tilt the 
world line, or the elastic string. 

2.2 Flux-line pinning in type-II superconductors 

The situation described by Eq. (5) is realized in type-II superconductors with 
extended defects. In a harmonic approximation, Nelson and Vinokur [29] de- 
rived a phenomenological Hamiltonian of a flux line in a superconductor with 
extended defects. The Hamiltonian takes the form of Eq. (5) with V(x) de- 
noting the pinning potential due to the defects. 

When an electric current is applied to a pure sample of a type-II superconduc- 
tor in the mixed phase, magnetic flux lines penetrating the sample are moved 
by electromagnetic forces, dissipate energy, and thus destroy the supercon- 
ductivity. Randomly located (but mutually parallel) extended defects such as 
columnar defects (typically created by bombardment of heavy ions) and twin 
boundaries (planer defects in the anisotropic YBCO) pin the flux lines effi- 
ciently as long as the flux lines are almost parallel to the defects [29-31]. When 
the external magnetic field is tilted away from the defects, but its transverse 
component H± is still small, we may expect that the bulk part of the flux line 
remains pinned (Fig. 2(a)). This is referred to as the transverse Meissner ef- 
fect [29], because the system in this region exhibits perfect bulk diamagnetism 
in the transverse direction. When one increases the tilt angle (Fig. 2(b)), a 
depinning transition occurs at a certain strength of the transverse magnetic 
field, namely H± c . (See Section 4 for further details.) 

The above depinning transition can be explained in a clear-cut way in terms 
of derealization in the random non-Hermitian quantum mechanics [1,2]. The 
depinning of the flux line corresponds to derealization of the relevant wave 
function of the quantum system; see Table 1 for other correspondence. 

Although it would be possible to describe the depinning within the framework 
of classical statistical mechanics, one can take advantage, in the non-Hermitian 
approach, of abundant results available concerning localization in Hermitian 
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(b) 



Fig. 2. Flux-line depinning due to tilt of the external magnetic field, (a) When the 
transverse component of the magnetic field is small, the flux line is pinned by a 
columnar defect in the bulk of the superconductor, although it is deflected near the 
surfaces, (b) For a larger H±, the flux line is depinned from the defect. 

Table 1 

Correspondence between the random non-Hermitian system and the flux-line sys- 
tem, 



Random non-Hermitian system 



Flux-line system 



(i-dimensional space 

and the imaginary time 
Quantum fluctuation 
World line 

Non-Hermitian field g 
Random potential V(x) 
Localization 
Derealization 
Real eigenvalue 

Complex eigenvalue in a periodic system 
Imaginary current of a quantum particle 



(d + l)-dimensional space 

Thermal fluctuation 
Flux line 

Transverse field H± 

Randomly located extended defects 

Pinning 

Depinning 

Transverse Meissner effect 
Helical structure of a flux line 
Tilt of a flux line 



random systems. It is also often practically easier to solve a Schrodinger equa- 
tion than to treat the same problem in the path-integral framework. 

There are many other physical realizations of the above random non-Hermitian 
system. Efetov [4] discussed the problem from the viewpoint of directed quan- 
tum chaos. Nelson and Shnerb found an interesting application to population 
biology [11]. In an independent work, Chen et al. [32] employed the same tech- 
nique to study sliding of charge-density waves in disordered systems. In the 
following, I concentrate on the flux-line analogy. 
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Fig. 3. The energy spectrum of the one-dimensional non-Hermitian model (2) with 
1000 sites. Each eigenstate e is marked by a tiny cross in the complex energy plane. 
Plots for different values of g are offset for clarity. The random potential at each site 
was chosen from a box distribution over the range [— t, t\. The same realization of the 
random potential {V x } was used for all plots here. A complex eigenvalue indicates 
that the eigenstate is delocalized and the corresponding flux line is depinned. 

3 Complex eigenvalues and derealization 

3.1 Energy spectrum in one dimension 

I now present an example of the energy spectrum of the random non-Hermitian 
system [1,2]. Figure 3 shows numerical results for the lattice Hamiltonian (2) 
in one dimension with 1000 sites. All energy eigenvalues are of course real for 
the Hermitian case g = 0. For weak g, e.g. g/Ti = 0.1 in the case of Fig. 3, 
all the eigenvalues are still real despite the fact that the Hamiltonian is non- 
Hermite. As we increase g, complex eigenvalues appear in the middle of the 
energy band and form a bubble. (The spectrum has an inversion symmetry 
with respect to the real energy axis, because the Hamiltonian (2) is a real 
matrix.) The region of complex eigenvalues expands towards the band edges 
as g is increased, and the whole spectrum eventually becomes almost elliptic. 
In fact, the spectrum exactly falls onto an ellipse if V x = [1,2]: 



Analytic forms of the random-averaged energy spectrum in one dimension have 
been obtained for weak disorder and weak g [7] and for general values of g with 
the Lorentzian random distribution [9,10]. By neglecting higher moments of 
the random distribution of the on-site potential {V x }, Brouwer et al. obtained 





an approximate shape of the bubble after averaging, in the form 

lime I = ^\/t 2 ~ (Ree) 2 , A " (7) 

& V 2^- (Ree) 2 



for | Ree| < e c and for small \g\, where A 2 = (V 2 ) is the second moment of 
the random distribution and 



The value of e c indicates the region of the bubble of complex eigenvalues and 
in fact is a "mobility edge" as discussed below. Note that the bubble does not 
exist for \g\ < hA 2 /(2t 2 ). 

For the Lorentzian distribution of the potential, P(V X ) = ^^^/(V 2 + 7 2 ), 
whose second moment is diverging, the exact shape of the bubble after aver- 
aging was obtained for general values of g [9,10]: 



Ree \ 2 /| Ime| + 7V 



cosh (g /ft) y \smh(g/h) J 



t 2 (9) 



for | Ree | < e c with the "mobility edge" 

e c = cosh(# /h)^Jt 2 - 7 2 sech 2 (#/ft). (10) 



The form (9) is remarkably simple; the upper and lower halves of the elliptic 
spectrum of the pure case, Eq. (6), are squeezed towards the real energy axis 
translationally by the distance 7, preserving the shape of the arcs. Eigenvalues 
that lost the support of the arcs become real and are distributed over the whole 
real axis except for the region of the bubble. The bubble entirely vanishes for 
\g\ < ftsinh~ 1 (7/t). 

Another approximate expression of the mobility edge e c has been obtained 
for general random distribution [6,13]. For the Lorentzian randomness, this 
expression is reduced to the exact result (10). 



3.2 Derealization and complex eigenvalues 



In the following, I explain the correspondence that is listed in Table 2. I 
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Table 2 

The delocalization criterion of the random non-Hermitian system. A wave function 
of the form ip ~ e - *^! for g = and \x\ — > oc gets delocalized when \g\ exceeds fin, 
and acquires a complex eigenvalue at the same time. 



non-Hermitian field g 


wave function 


eigenvalue 


\g\ < Hk 


localized 


real 


\g\ > Tin 


delocalized 


complex 




Fig. 4. The world line (the thick line) of a current-carrying particle forms a helix 
when periodic boundary conditions are imposed in the space direction. Hence peri- 
odicity appears in the imaginary-time direction. The periodicity is described by an 
oscillatory factor e lrIme . 

first argue that a complex eigenvalue indicates a delocalized wave function. 
Consider the imaginary-time dynamics of a wave function 

ip(x; t) = ip(x) e~ T£ oc e - irIm£ , (11) 



where e is the eigenvalue of the eigenfunction i]){x). The above equation shows 
that the imaginary part of the eigenvalue, if it exists, gives rise to an oscillatory 
behavior of the system in the imaginary-time direction. In fact, the oscillatory 
behavior comes from the world line wrapping around the system as a helix 
(Fig. 4). Note here that periodic boundary conditions in the space directions 
are imposed on the quantum system. Hence a quantum particle, when it is 
delocalized, circulates around the system. This yields a helical flow ascending 
in the imaginary-time direction of the (d + l)-dimensional imaginary-time- 
space. The pitch of the helix corresponds to the reciprocal of the imaginary 
part of the energy. 
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It follows from the above argument that a wave function of the periodic system 
acquires a complex eigenvalue as soon as it gets delocalized, or the correspond- 
ing flux line gets depinned and tilted. Thus observing the energy spectrum of 
the non-Hermite Hamiltonian is a convenient way of investigating the flux-line 
depinning. 

Next, I explain the derealization criterion \g\ 5 Tin in Table 2. For this pur- 
pose, I introduce the imaginary gauge transformation [33]. Suppose that we 
obtain an eigenfunction of the Hamiltonian for g = with a real energy eigen- 
value: 

Ho^o(x) = € ipo(x). (12) 

Since the non-Hermitian field g is equivalent to an imaginary vector potential, 
we may be able to gauge out the field from the non-Hermite Hamiltonian H.(g). 
In other words, the eigenfunction of the Hamiltonian H(g) corresponding to 
■00 m &y be given by 

iP(x) = e^ /n ip (x) (13) 
and the eigenvalue eo may remain the same. 

The above transformation is valid only in a certain range of g. Assume that 
the wave function for g = is asymptotically given by i^o(x) ~ e -re l x l, where 
k is a constant. Equation (13) then takes the form 

ip(x) ~ e - K W + * s/h . (14) 

This wave function is (asymmetrically) localized for \g\ < Tin. In this region 
the function (14) satisfies periodic boundary conditions asymptotically in the 
infinite-system-size limit. Hence the imaginary gauge transformation is valid 
for \g\ < fiK and Eq. (13) is indeed the eigenfunction of the non-Hermite 
Hamiltonian 7i{g) with the real eigenvalue e . In fact, closer inspection of the 
spectrum in Fig. 3 would reveal that the real eigenvalues in the "wings" of the 
spectrum does not depend on g. This rigidity reflects the transverse Meissner 
effect of the corresponding flux-line system. 

For \g\ > Tik, on the other hand, the function (14) does not satisfy periodic 
boundary conditions, because it blows up in the direction of g in this region. 
Thus Eq. (13) is no longer an eigenfunction of the non-Hermite Hamiltonian 
Ti.(g). It is in this region that the eigenfunction is delocalized and acquires 
a complex eigenvalue. (Equation (13) is an eigenfunction if we impose open 
boundary conditions. All the eigenvalues remain real in this case, although 
the eigenfunctions are nonetheless delocalized. This is consistent with the fact 
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that the helix in Fig. 4 never appears if the boundaries are open in the spatial 
directions.) 



3.3 Mobility edge and inverse localization length 



According to the above argument, the eigenstates that belong to the bubble 
in Fig. 3 are delocalized, while those in the wings of the real eigenvalues are 
localized. Hence the two vertices of the bubble, e c and mobility edges. 

The complex energy eigenvalues appear first in the middle of the energy band 
because, in the one-dimensional Anderson model (g — 0), k is the smallest for 
the eigenstate at e — 0, as exemplified below. 

Table 2 provides a convenient method of estimating the inverse localization 
length for g — 0. It is deduced from the derealization criterion that a state 
at a mobility edge for a value go of the non-Hermitian field has the inverse 
localization length k — \go\/h for g = 0. A numerical calculation of k by this 
method was presented in Ref. [2] for the one-dimensional lattice model with 
a box distribution. Using the analytic result (10), we can also calculate the 
inverse localization length of the Lloyd model [34] (the one-dimensional An- 
derson model (g = 0) with the Lorentzian random distribution) as a function 
of the energy, by solving 

2 ' 7 \ 2 +2 



e = cosh/tyt 2 - 7 2 sech 2 /t, or — - — + — — = t. (15) 

\cosh/c/ Vsinh/t/ 



The solution 



,( ()=aA -.MiMl^ (16) 



reproduces the exact result obtained by Hirota and Ishii [35,36] and Thou- 
less [37]. 

Figure 5 shows the function k{s) for t — 7 = 1. If we apply to the system the 
non-Hermitian field of the value, say, g/h = 1.5 as indicated by a dotted line 
in Fig. 5, the eigenstates below the dotted line (— e c < e < e c ) get delocalized 
and form a bubble in the spectrum, while those above the line (e < —e c and 
e > e c ) remain localized and stay in the wings of the spectrum. 
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Fig. 5. The inverse localization length (16) of the Lloyd model (the solid line). The 
parameter values used in this figure are t = 7 = 1. The crossing points of the solid 
line and the dotted line (indicating a value of g/fi) yield the mobility edges e c and 

Sc. 

3.4 Other one- dimensional models 

Feinberg and Zee [8,12] introduced an interesting limit of the lattice Hamilto- 
nian (2) in one dimension, namely the one-way model: 

H = Y,(-\x + l){x\+V x \x)(x\). (17) 

x 

In other words, they took an extremely non-Hermitian limit, te 9 ^ n — > 2 and 
te~ 9 / n — > 0. Exact calculations of the spectral curve (Regime) and the mo- 
bility edge e c become possible for various random distributions including the 
box distribution, the binary distribution and even diluted randomness. See 
Ref. [12] for details. 

In Figure 6, I show a numerical result of the energy spectrum of a ladder 
system, 

H = H 1 + H 2 -lj2(\x,2)(x,l\ + \x, 1) (x, 2|) , (18) 



where Tii and is the non-Hermitian Anderson Hamiltonian for each leg of 
the ladder while the last term denotes hopping between the legs; "x, 1" and 
u x, 2" denote a site on the first and second leg, respectively. Alternatively we 
may regard the labels 1 and 2 as an additional degree of freedom such as spin 
or flavor. On the basis of the argument in the previous subsection, we can 
deduce from the result in Fig. 6 that there are two minima of k(e) in the 
Hermitian ladder system. 



12 




Fig. 6. The energy spectrum of the non-Hermitian ladder model (18) with 500x2 
sites. Each eigenstate e is marked by a tiny cross in the complex energy plane. The 
random potential at each site was chosen from a box distribution over the range 
[— t,t]. The same realization of the random potential {V x } was used for all plots 
here: (a) g/K = 0.3, (b) g/h = 0.5, (c) g/h = 0.7 and (d) g/h = 0.9. 

The above numerical calculation was immediately followed by Zee's analytic 
calculation of the spectrum for the Lorentzian randomness [14]. Just as in the 
one-dimensional case (10), the averaged spectrum is squeezed towards the real 
axis as the randomness is increased. Thereby we can exactly calculate k(e) for 
the ladder Lloyd model. The result is simply superposition of two functions of 
the form (16). This is consistent with the above numerical result for the box 
distribution. 

Note that, except in the case of the single chain, k(e) thus calculated is an 
upper bound of what is referred to as the "inverse localization length" in 
the context of the Anderson localization [38-40]. The geometric average of 
the Green's function is taken in denning the Anderson-localization length, 
while the arithmetic average is taken in the above analytic calculation. (The 
argument that yielded Table 2 is still valid for each realization of the random 
potential.) 
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3.5 Energy spectrum in two dimension 

Finite-size calculations of the energy spectrum of the non-Hermite Hamilto- 
nian (2) for d = 2 [1,2] appeared to suggest the following three regions of 
the non- Hermit ian field g. First, it is widely accepted for g = that all eigen- 
states of two-dimensional random systems are localized with finite localization 
lengths. Hence, by using the imaginary gauge transformation again, we can 
conclude that there is a finite region of small g where all states remain lo- 
calized. This is consistent with the finite-size data [1,2]. As g is increased, 
delocalized states with complex eigenvalues appear as in the one-dimensional 
case. For intermediate values of g, however, the energy spectrum shows much 
more complicated structure than in the one-dimensional case; see Ref. [2] for 
details. For larger g, the spectrum becomes similar to the one of the system 
without impurities, as was the case in d = 1. 

Nelson and Shnerb [11] suggested that the third region of g disappears in the 
thermodynamic limit. For a very large g, or a very large transverse magnetic 

— * 

field H±, the flux line lies nearly sideways in the superconductor. When we 
project the flux line and the columnar defects onto a two-dimensional plane 
from above, the problem is approximately reduced to a string lying on a plane 
with point impurities. Nelson and Shnerb analyzed this reduced problem in 
terms of Burger's equation with noise [41] (or the Kardar-Parisi-Zhang equa- 
tion [42]) and predicted a fractal geometry of the flux-line configuration. This 
fractality may cause the complicated energy spectrum observed in the sec- 
ond region of g. The fractal geometry does not emerge until the system size 
becomes larger than a certain crossover length, which can explain the appear- 
ance of the third region of g in the finite-size data [1,2]. The above theoretical 
prediction is consistent with numerical results in Ref. [32] , although the numer- 
ical estimate of the exponent characterizing the fractal geometry is somewhat 
different from the theoretical value predicted by Nelson and Shnerb [11]. 

Note, however, that a different conclusion may be suggested by Zee's ana- 
lytic calculation of the two-dimensional spectrum for the Lorentzian random 
potential [14]. In the case of the Lorentzian random distribution, just as in 
one dimension, the averaged spectrum is squeezed towards the real 
the randomness is increased. In the process, the spectrum keeps the regular 
structure of the spectrum of the non-random system. 



4 Transverse Meissner effect 

One of the interesting issues in the context of the flux-line depinning is how the 
transverse Meissner effect breaks down as the depinning point is approached 
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from the side of the pinning phase. As is shown schematically in Fig. 2, the 
pinning breaks first near the surface. This is also observed in numerical re- 
sults [2]. Thus we can define a penetration depth r* for the transverse Meissner 
effect as the typical thickness of the sub-surface region where the flux line is 
deflected from the pinning center. (Note that r* is different from the pene- 
tration depth of the underlying superconductivity.) The penetration depth r* 

— * 

diverges with a certain exponent, as we increase H±, or g. This divergence 
leads to the breakdown of the bulk pinning. The conclusion of Refs. [1,2] is 

T*~(H Lc -H ± )~ d , (19) 

where H± c is the depinning field. Note that the dimensionality of the super- 
conductor is d + 1. 

Since we regard the flux line as the world line of a quantum particle in the 
framework of the non-Hermitian quantum mechanics, the deflection of the flux 
line is interpreted as hopping of the quantum particle from a strong impurity 
(the localization center) to weaker impurities. Near the derealization point 
g c (oc H± c ), Mott's argument of variable-range hopping [43] is readily appli- 
cable even to the non-Hermitian case. Minimizing a hopping matrix element 
near the surface leads to Eq. (19). 



5 Summary 

A class of random non-Hermitian quantum systems has been found to be 
relevant to various physical systems and to show intriguing properties. In par- 
ticular, the derealization phenomenon of the random non-Hermitian system 
is equivalent to the depinning of a flux line in type-II superconductors. We 
can investigate the derealization just by observing the energy spectrum of 
the non-Hermite Hamiltonian. Various depinning phenomena may be under- 
stood in a similar way. Interesting future problems include many-band non- 
Hermitian Anderson models, detection of the fractality in the spectrum of 
the two-dimensional system, and generalization of the theory to the case of 
interacting systems. 
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A Path-integral mapping 



I here show derivation of Eq. (4) from Eq. (3). There are various ways of for- 
mulating the path-integral mapping. The following is based on the formulation 
in Ref. [28], which employs the Trotter decomposition [15]: 

e-W* = hm ( e -VAr/n e -KAr/ny^ (A l) 

Here /C and V denote the kinetic and potential terms of the Hamiltonian (1), 
respectively, and At = L T /n. 

The matrix element (3) is transformed as follows: 



2 = Jim <y | ( e - VAT / R e- /CAT / R ) n | ^ 

- n - n 



k=0 



k=l 



x n i ex p 

fc=i l 

The plane wave is given by 



x, 



) exp 



At 
2mh 



Pit 



(27r^) d / 2 



exp 



TVk ■ Xk 
n 



Xk-l 

(A.2) 
(A.3) 



Hence the integral over each pk in Eq. (A.2) is a Gaussian integral, 



dp k exp 



^-(p k + ig) 2 + ^P k -Ax k 



m \ d l 2 
:) ex P 



m 



—q- Ax k — 

7i 2KAt 



(Ax k ) 2 



^hAT, 

where Axk = %k — %k-i- Thus Eq. (A.2) is reduced to 



(A.4) 



/ 777 \ nd/2 r ( n \ 



X exp 



At " 



m ( Ax*k 
~2 1^7 



_ Af fc 

5 ' a7 + v{xk) 
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m / dx(r) 
2~ [ dr 



dx(r 
dr 



+ V(x(r)) 



(A.5) 



By assuming the free boundary conditions ^'(x) = ip (x) 
at Eq. (4) with Eq. (5). 



const., we arrive 
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